First Order Superfluid to Bose Metal Transition in Systems with Resonant Pairing 
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Systems showing resonant superfluidity, driven by an exchange coupling of strength g between 
uncorrelated pairs of itinerant fermions and tightly bound ones, undergo a first order phase transition 
as g increases beyond some critical value gc- The superfiuid phase for g < (?c is characterized by 
a gap in the fermionic single particle spectrum and an acoustic sound-wave like collective mode of 
the bosonic resonating fermion pairs inside this gap. For g > gc this state gives way to a phase 
uncorrelated bosonic liquid with a spectrum. 
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Introduction. The issue of the crossover between a BCS 
type superfluid and a Bose Einstein condensation (BEC) 
of tightly bound fermion pairs has received increasing 
attention ever since it became feasible to experimentally 
control such a crossover in fermionic atomic gases in op- 
tical traps and lattices[l]. From the theoretical side, such 
a crossover had until then been mainly studied for single 
component fermionic systems with static, short-range at- 
tractive inter-particle interaction[2]. It became clear that 
in these systems one is confronted with a unitarity limit 
problem, where the scattering length abruptly changes 
from positive to negative infinity. Numerical as well as 
analytical results [3] suggest that around such a unitarity 
point, one should expect, effectively, two kinds of quasi- 
particles coexisting with each other: itinerant fermions 
and tightly bound fermion pairs behaving as bosons. 

Such a situation had been anticipated in connexion 
with the physics in electron-lattice coupled systems in a 
regime between weak and strong coupling and resulted 
in the proposition of the boson- fermion model (BFM) [4] . 
The essence of the BFM scenario is to assume itin- 
erant fermions and localized bosons (made up of two 
tightly bound fermions) - the two species which repre- 
sent the physics of the extreme weak and strong cou- 
pling limits of such systems. The crossover behavior be- 
tween those two limiting regimes is then parametrized 
by an exchange coupling of strength g between the two 
species. This scenario can be justified for: (i) fermionic 
atomic gases, where pairs of fermions exist in different 
electron- nuclear spin configurations, favoring or disfavor- 
ing binding between them (magnetic field tunable Fes- 
hbach resonance) [5] , (ii) hole pairing in strongly cor- 
related systems, showing plaquette RVB states[6] and 
(iii) resonating bipolarons in moderately strongly cou- 
pled electron-lattice systems [7], for which this scenario 
was anticipated in the first place [4]. 

Contrary to fermionic systems with attractive inter- 
particle interaction, where the ground state is always 
superfluid, the BFM scenario shows a superfluid (am- 
plitude driven) BCS-like state for weak coupling, while 
for strong coupling it represents a liquid of spatially 



phase uncorrelated bonding pairs: 1 / A/2(uiC-^c^ -I- w,;6j). 

= \u,\e^*^/^v, = \v,\e-"f>^/^ (with \u,\^ + \v,\^ - 1) 
denote the local amplitudes and phases of the two com- 
ponents, and c|g., respectively b], the creation operators 
of the itinerant fermions (with spin a) and of the local- 
ized bosonic tightly bound fermion pairs at site i. The 
emergence of a macroscopic superfluid state requires the 
onset of long-range phaselocking of such bonding pairs. 
This needs a concomitant weakening of the local phase 
correlation of the bonding pairs [8]. It is this competition 
between local and global phaselocking which is at the ori- 
gin of a phase transition between a superfluid state for 
<? < ffc and a state of spatially phase uncorrelated bond- 
ing pairs for g > gc- In this Letter, we argue that this 
transition is of first order. 

The exact nature of that phase uncorrelated state is 
presently not yet fully understood. For very large val- 
ues of g, it eventually must be an insulator. Here, we 
show that close to gc it constitutes a Bose Metal made 
out of such bonding pairs and present the first study 
of the excitation spectra expected for such a BFM sce- 
nario, which treats the fermions and bosons on equal 
footing as we vary the exchange coupling g. This will 
be done on the basis of a renormalization procedure for 
IIamiltonians[9, 10] designed in such a way as to op- 
timally provide a fixed point Hamiltonian with single- 
particle spectra for effective fermionic and bosonic quasi- 
particles and reducing the interactions between them to 
higher order correction terms such that, in a first step, 
they can be neglected. 

The model. The Hamiltonian for the BFM is divided 
into a free and a boson-fermion exchange coupling inter- 
acting part, i.e., H — Hq + Hint- 

Ho = J2^4- + - 2^)6^6, 

k.cr q 

Hint = ^i9k,pbl^pCk,lCp,^ + gl^pbk+pcl -fcl^^) (1) 

Important for this model is to remember that the 
bosons are composed of two fermions. This is taken care 
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FIG. 1: (Color online) Phase diagram of the BFM. The thin 
solid line marks the transition from a Fermi liquid to a BCS- 
like superfluid phase which merges into a crossover (shaded) 
with the pseudogap phase. The thick solid (dashed) line mark 
the first (second) order transition from a superfluid of locally 
phase correlated boson-fermion bonding pairs to an itinerant 
bosonic metal (Bonding Pair Bose Metal). Diamonds repre- 
sent the results of the present calculation. 



of by requiring that both types of particles have the same 
chemical potential ^, fixed such that the total particle 
density ntot = npi + npi + "^ns is constant. We shall in 
the present study consider the case of low overall density 
and choose for that purpose nF,o- — ns — 0.25. 

Fig. 1 shows the generic phase diagram of the BFM at 
low densities. For small values of exchange couplings g, 
the thin solid line indicates the transition from a Fermi 
liquid (FL) to a BCS type superfluid phase, controlled by 
amplitude fluctuations and characterized by the opening 
of a gap in the single particle fermion spectrum. This 
transition line merges, with increasing g, into the pseu- 
dogap crossover (shaded) which separates a canonical FL 
from a fermionic system with strong pair correlations 
(pseudogap phase). This crossover was analyzed in de- 
tail in Refs. 11, 12. In the present work, we focus on the 
large exchange coupling regime and the transition from 
a superfluid to a Bose Metal phase of resonating fermion 
pairs. 

The method. The flow equation approach consists in 
performing continuous unitary transformations upon the 
initial Hamiltonian until the Hamiltonian is (block) di- 
agonal. This is done in a systematic way by first decou- 
pling states with high energy differences and, in the end, 
decoupling almost degenerate states. Previous applica- 
tions of continuous unitary transformations considered 
the weak coupling regime of the BFM, where the zero 
momentum Cooper pair channel is the dominant interac- 
tion and the bosonic dispersion can be approximated by 
its q = 0-value[ll]. In order to describe the evolution of 
the Hamiltonian in the strong-coupling phase with and 



without global phase coherence, it is crucial to treat the 
bosonic dispersion in the whole Brillouin zone. 

Flow equations are formulated in terms of an anti- 
Hermitean generator 77(^) and read 9£if(£) = [r]{£), H{£)], 
with i being the flow parameter. In this work, we 
shall choose the generator canonically, i. e., ri{£) = 
[Hoie),H{e)], implying a,Tr(iJ(£) - iJoW)' < and 
hence gk,p — > for — > oo [10]. The infinitesimal trans- 
formations induce a flow on the fermion as well as boson 
dispersion and exchange coupling constants gk,p- 

degk^p = -al,p9Kp, ak^p = 4 + 4 ^ ^'=+p 

= ^ 51 ("p,fel5'p,fcP'5^,T + ak,p\gk,p?Sa^i) 
p 

diEq = — ^{aq^pjgq^pj'^[-l + n[^fj^] 
p 

(^p.,q-p\9p.,q-p\ '^\^qJp\ (2) 

For £ — oo, the renormalized fermions and bosons are de- 
coupled and characterized by the fixed point dispersions 
e* = efc(£ = oo) and E* = Eq{i = oo). 

In order to close the infinite hierarchy of newly gener- 
ated interaction terms, higher order interaction terms in 
their normal ordered form with respect to the fixed point 
Hamiltonian Hq are neglected. To lowest order, this in- 
troduces bilinear expectation values of the fermionic and 
bosonic operators. We will choose the expectation values 
explicitly ^-dependent, i. e., 

nf ^)(^) = {b%){e) = (e(i^,W-2M£))/fc«T -iy\ (3) 

The chemical potential fi{£) is determined such that 
the total number of particles iVtot = J2ka ^io-^^ (^) + 
rig {t) is conserved, which induces a flow of ^{i) 
with a fixed point value ^* = ^{l = oo). 

The justification for assuming fermion and boson dis- 
tribution functions, Eq. (3), in form of the distribution 
functions for non-interacting quasi-particles is the follow- 
ing: For small £, states with high energy differences are 
decoupled, which corresponds to the perturbative regime 
and thus to well defined fermionic and bosonic quasi- 
particles. For large the exchange interaction g^^p has 
basically dropped to zero, according to how this renor- 
malization procedure was constructed. It is this asymp- 
totic regime £ — > cx) which determines the critical be- 
havior of the system [13]. Since the main effect of the 
exchange coupling g is not to induce a significant change 
in the relative number of fermions or bosons but to in- 
troduce phase coherence between the Cooper pairs and 
bosons, Eq. (3) represents a good approximation also in 
the strong coupling regime [14]. Hence, even though the 
truncation scheme is perturbative in the exchange cou- 
pling, the strong coupling regime can be discussed within 
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the flow equation approach (as has been shown also for 
the Hubbard model, see e.g. Ref. 16). 

Numerical results. As initial conditions for the renor- 
malization flow, we choose a dispersionless bosonic band 
with Eq = Ab = —0.6 and a fermionic tight-binding 
dispersion Ck.a = Cfc = —2tcoa{ka) where we set the lat- 
tice constant a = 1 and use the bandwidth D — it a.s 
energy unit. We further assume a local exchange inter- 
action gk^p = g. For the numerical integration of the flow 
equations, Eq. (2), we choose a one-dimensional system 
with N = 100 lattice sites which already resembles the 
thermodynamic limit. [15] As we shall discuss below, the 
phase transition is determined by the chemical potential 
moving out of the renormalized fermionic band of the 
fixed point Hamiltonian. Since the divergencies of the 
density-of-states at the bare band edges is smeared out 
in the course of the renormalization flow, the basic re- 
sults are expected to remain qualitatively valid also in 
higher dimensions. 

The present study, concentrating on the strong cou- 
pling regime, deals with the nature of transition from the 
superfluid state of resonating fermion pairs to their phase 
uncorrelated Bose Metal state. The critical interaction 
strength gc as function of temperature T is indicated by 
diamonds (<)) in Fig. 1. Such a Bose Metal consists of 
current carrying states being composed of bonding pairs 
without any contributions from single particle fermion 
states. The transition is characterized by fi* moving 
out of the fermionic band e^. It is discontinuous for 
T < Tq ~ 0.75 (thick solid line in Fig. 1) and changes 
into a continuous transition for T > Tq (dashed line). 
The nature of the transition can be deduced from the 
renormalized fermionic dispersion which is shown in 
the upper panels of Fig. 2 for various exchange interac- 
tions g and temperatures T — 0.01 and T = 0.1. 

For T = 0.01 and g < g^ ~ 0.8 (sohd Unes of the 
upper left panel of Fig. 2), the fermionic excitations 
show a single band separated by a gap which is rem- 
iniscent of a lower and upper Bogoliubov band, when 
keeping only the parts with maximal spectral weight. As 
g increases, the gap gets bigger and the fraction of co- 
herent fermions which are separated from incoherent ex- 
citations by the superconducting gap becomes smaller. 
For g > gc (dashed lines), the lower Bogoliubov band 
has disappeared and the upper Bogoliubov band now 
plays the role of our effective fermionic band which 
has abruptly moved up in energy and now lies above ^* . 
The fermionic spectral function now consists of purely in- 
coherent contributions characterized by breaking up the 
strongly bound fermion pairs in the purely local-phase 
correlated liquid. The total number of coherent and in- 
coherent fermions is approximately constant since the ex- 
change coupling g does not significantly change the rel- 
ative number of bosons and fermions. For T = 0.1, the 
change of the fermionic dispersion is continuous for in- 
creasing g when fi* comes to lie below e^^g, as shown on 
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FIG. 2: (Color online) The fermionic and bosonic fixed 
point dispersion (upper panels) and Eg (lower panels). 
Left panels: At T = 0.01 for g = 0.05 (solid bold line), 
g = 0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.6,0.7,0.79 
(solid line) and g — 0.8,0.9 (dashed line). Right panels; At 
T = 0.1, for g = 0.25,0.32,0.33,0.34,0.35 (solid line) and 
g = 0.36, 0.37, 0.38, 0.39, 0.4, 0.45 (dashed line). 



the upper right hand side of Fig. 2. This yields gc — 0.35. 

At low temperatures and close to the phase transition, 
we notice a qualitative change in ^* as a function of g for 
different temperatures. Fig. 3 shows ^* as a monoton- 
ically decreasing function of g for T = 0.1, changing to 
a non-monotonic behavior at T = 0.5. This is indicative 
for the continuous transition for T > Tq ~ 0.75 changing 
into a first order transition for T < Tq where we find the 
coexistence of two phases with different relative densities 
of bosons and fermions. This feature is accompanied by 
a smooth, respectively abrupt, change of the position of 
fi* from inside the fermionic band for g < 5c to outside 
of that band for g > gc, as seen in the insert of Fig. 3, 
where e^^q — /i* is plotted as function of g. 

For T < Tq, also the bosonic dispersion changes 
abruptly from a linear to a quadratic spectrum at low 
q as g increases beyond gc (lower left panel of Fig. 2). 
The initially dispersionless bosons thus first become su- 
perfluid and then itinerant as g increases from weak to 
strong coupling. For Tq < T — 0.1 (lower right panel), 
the transition is continuous, nevertheless the width of 
the occupation number n^^(£ = oo) around q — oi the 
renormalized dispersion changes rapidly, indicating the 
transition from superfluidity to a Bose metal at finite 
temperatures. 

Summary and discussion. We have discussed the 
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FIG. 3: (Color online) The renormalized chemical potential 
fi* as function of the coupling strength g with respect to the 
critical values and gc of the phase transition from superflu- 
idity to a Bose Metal. The curves correspond to three different 
temperatures T > To, T « To and T < To where To defines 
the transition from second to first-order phase transition. 



boson-fermion model within the flow equation approach, 
including both the Cooper-pair and the density-density 
channel for all wave numbers. We find a first order 
phase transition from a superfiuid phase to a Bose Metal 
for low temperatures in the strong exchange coupling 
regime. In the fermionic one-particle spectral function, 
both phases show a gap. But, depending on the rela- 
tive strength of the hopping matrix element t with re- 
spect to the g, this leads either to a phase-correlated or 
to a phase-uncorrelated bonding-pair liquid. The transi- 
tion is seen in the bosonic dispersion as a change from 
linear to quadratic in momentum behavior, i.e., in the 
superconducting regime we have the Goldstone mode 
E* — E*^Q = vsfQ (for small q), whereas in the phase- 
uncorrelated regime we have E* — E*^q = (f' jlm* with 

In the superfiuid phase, the renormalized dispersion for 
the bosons is linear for wave numbers g = up to some 
9max which depends on the size of the regions where local 
phase correlations between the itinerant fermions and the 
tightly bound boson pairs prevail. They become stronger 
for g ^ gc which renders impossible the maintenance of 
the spatial phase correlations for the Cooper-pairs in the 
fermionic subsystem over short distances and hence a lin- 
ear dispersion of the bosons for large momenta. Qmax thus 
decreases with increasing g and eventually reaches zero 
at the phase transition. We note that the slope of the lin- 
ear dispersion, i. e., the superfiuid velocity wsf, hardly 
changes even for g gc- This is yet a further manifesta- 
tion of the first order phase transition (in case of a second 
order phase transition the slope would eventually tend to 
zero) . 

In the strong-coupling Bose metal phase [g > gc), all 
one-particle fermion states exclusively exist as compo- 



nents of the bonding pairs. The single-fermion particle 
spectrum then exhibits a correlation gap and the strong 
local phase coherence, leading to a locally fiuctuating 
field, destroys all coherence of the itinerant fermions. 
The bosonic dispersion, E*, of the renormalized Hamil- 
tonian corresponds, within such a fiow equation pro- 
cedure, to renormalized bosonic operators of the form 
Uqbl + Efc^fc.gcI+g.T^tfc.i [12]. Since E* oc q^, it rep- 
resents a Bose Metal composed of the above mentioned 
bonding pairs. 

The first order transition has been deduced from the 
dynamical properties of the system. In order to discuss 
the transition within thermodynamic quantities, one has 
to study the fiow of the operators in addition to the one 
of the Hamiltonian. This is a separate problem which will 
be addressed in the future. The present study dealt with 
spatially homogeneous systems, but should have its bear- 
ing on inhomogeneous systems, such as fermionic gases 
in optical traps. 
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